Numerical assessment of irreversibility in radiated Sutterby nanofluid flow with activation energy and Darcy Forchheimer

Entropy generation is a concept that is primarily associated with thermodynamics and engineering, and it plays a crucial role in understanding and optimizing various processes and systems. Applications of entropy generation can be seen in turbo machinery, reactors, chillers, desert coolers, vehicle engines, air conditioners, heat transfer devices and combustion. Due to industrial applications entropy generation has gained attention of researchers. Owing such applications, current communication aims to model and analyzed the irreversibility in Sutterby nanoliquid flow by stretched cylinder. Momentum equation is reported by considering porosity, Darcy Forchheimer and magnetic field. While in energy equation radiation and Joule heating effects are accounted. Activation energy impact is accounted in the modeling of concentration equation. Thermodynamics second law is utilized for physical description of irreversibility analysis. Through similarity transformations dimensional equations representing flow are transformed to dimensionless ones. Numerical solution for ordinary system is obtained via Runge–Kutta-Fehlberg scheme in Mathematica platform through NDsolve code. Influence of prominent variables on velocity, entropy, temperature, Bejan number and concentration are graphically analyzed. Coefficient of skin friction, gradient of temperature and Sherwood number are numerically analyzed. The obtained results show that velocity field decreases through higher porosity and Forchheimer variables. Velocity and temperature curves shows an opposite trend versus magnetic parameter. A decay in concentration distribution is noticed through larger Schmidt number. Entropy generation amplifies against magnetic parameter and Brinkman number.


List of symbols m
Power law index (−) k Thermal conductivity W m −1 K −1 τ Heat capacity ratio (−) B 0 Strength of magnetic field kg s Density of fluid kg m −3 (u, w) Velocity components ms −1 T Temperature (K) Brownian motion coefficient m Viscoelastic liquids like molten polymers and polymer solutions exhibit numerous verities of rheological characteristics.These features include time dependent viscosity, shear rate dependent viscosity, normal stresses in steady shear flow and various time-dependent elastic effects.Sutterby fluid model is one of the most important models proposed by Sutterby 1 , which addresses the solution of high polymer aqueous solutions.Imran et al. 2 inspected the thermal radiation and chemical reaction influences on Sutterby nanomaterial flow in inclined elastic channel.Impact of radiation, mixed convection and Joule heating on Sutterby nanoliquid flow in a vertical channel is numerically examined by Hayat et al. 3 .Akbar and Nadeem 4 reported mixed convection flow of Sutterby nanoliquid in a diverging tube.Helical flow of Sutterby nanomaterial between two concentric cylinders is scrutinized by Batra and Eissa 5 .Ishtiaq et al. 6 examined shear thickening/thinning behavior of Sutterby nanomaterial flow over biaxially stretchable sheet with magnetic field and heat source/sink.Mixed convection and Arrhenius kinetics effects on 3-D steady flow of Sutterby nanomaterial is conveyed by Azam et al. 7 .Khan et al. 8 inspected the features of stratified Sutterby nanomaterial flow in presence of external radiation and Lorentz force.Mixture of nano-sized metallic particles and base fluids are nanofluids/nanomaterials.Nano-sized metallic particles include metallic oxides, metals, carbon nanotubes and nitrides.Conventional base fluids are water, ethylene glycol and light oils.Nanofluids have higher thermal performance as compared to conventional carrier liquids.Practical usages of such nanofluids can be seen in thermal engineering processes like fuel cells, refrigerators and engine oil.Firstly Choi 9 added nano-sized metallic particles in carrier fluids and concluded that thermal features improved significantly.Buongiorno 10 provided a model to study heat transfer augmentation in nanomaterials.He considered seven slip mechanisms for nanoparticles and proved that thermophoresis diffusion and Brownian movement are governing factors as compared to others.Prasad et al. 11 reported radiative nanofluid flow with Lorentz force effect.Turkyilmazoglu 12 inspected mass and heat transportation in nanofluid flow over different frames using Buongiorno model.Tian et al. 13 analyzed convectively heated MHD nanoliquid flow having stagnation point over stretching sheet.Features of laminar flow of viscous fluid flow by stretchable cylinder are analytically examined by Turkyilmazoglu 14 .Hayat et al. 15 analyzed convective hybrid nanofluid flow with radiation and heat transfer characteristics.Influences of Hall current and electrical MHD in flow of micropolar nanofluid between a pair of rotating disks is explored by Awan et al. 16 .Hussain et al. 17 inspected features of rotating flow of hybrid nanoliquid accounting the influences of restricted slip boundary constraints.Qureshi et al. 18 examined impressions of heat generation and magnetic field in hybrid peristaltic flow in a metachronal wave.Parveen et al. 19 inspected the characteristics of dissipative bioconvective flow of nanoliquid which contains chemotactic microorganisms through Joule heating.Rheological properties of Pseudo plastic nanomaterial flow in a symmetric channel accounting the effects of ciliary motion is presented by Khan et al. 20 .
Bayesian regularization networks is implemented to explore the features of nanofluid flows considering various effects by Awan et al. 21,22 .Ahmad et al. 23 deliberated Maxwell nanomaterial flow by permeable rotating frame with heat transfer analysis.Flow behavior of radiative MHD Maxwell nanoliquid considering viscous dissipation effects is investigated by Hsiao 24 .Hassan et al. 25 inspected the behavior of MHD hybrid nanomaterial flow which contains SWCNT-Ag as nanoparticles with radiation effects.7][28][29][30] .
Due to applications in environmental, chemical, industrial and pharmaceutical sectors nanoliquid flows over porous surface have gained significance.Applications of such flows may include in energy storage units, geothermal heat exchanger layouts, nuclear waste disposal and crude oil production.Seddeek 31 reported convective nanofluid flow having thermophoresis and dissipation features saturating porous space.Umavathi et al. 32 described Darcy-Forchheimer convective nanofluid flow with Brinkman relation.Muhammad et al. 33 analyzed Maxwell nanoliquid flow in Dacry-Forchheimer porous space.Hayat et al. 34 scrutinized Darcy-Forchheimer flow with carbon nanotubes in presence of permeability.Alzahrani 35 investigated bidirectional flow of CNTs in Darcy-Forchheimer and porous surface.Turkyilmazoglu 36 studied the features magnetic field of uniform strength applied horizontally to the flow generated by rotating disk.Hayat et al. 37 discovered Carreau nanoliquid flow due to heated surface with Darcy-Forchheimer and porosity effects.
In chemical reactions role of activation energy (AE) is noteworthy.The lowest energy amount obligatory to trigger a chemical reaction is named as AE.In different fields AE concept is widely used such as water mechanics, oil emulsions and oil tank counting.Bestman 38 revealed flow over porous surface in existence of energy activation and binary response with mass and heat transmission characteristics.Makinde et al. 39 reported steady radiative flow over porous medium with heat transfer and chemical reaction impacts in an optically reedy atmosphere.Maleque 40 discovered unsteady boundary layer flow having heat sink/source and energy activation impacts.Awad et al. 41 investigated unsteady rotational viscid fluid flow considering AE and chemical reaction.Impacts of AE and thermal radiation on MHD flow of Carreau are reported by Hsiao 42 .
Entropy generation (EG) is an extensive property of thermodynamics.Thermodynamics second law states that in an isolated thermal system entropy never decays.In an irreversible reaction total entropy always increases while in reversible reactions it remains constant.Bejan 43 studied irreversibility in convective nanoliquid flow.Turkyilmazoglu 44 scrutinized EG and slip impressions in radiated fluid flow in by metallic permeable channel.Khan et al. 45 calculated the EG in radiative flow of Sisko liquid by stretched surface with dissipation effect.Vatanmakan et al. 46 examined EG in steam flow with volumetric heating and turbine blades.Hayat et al. 47 inspected EG in Casson type nanoliquid flow by stretchable sheet through magnetic field and Arrhenius kinetics.Gul et al. 48eported EG in viscoelastic Poiseuille nanofluid flow.Xie and Jian 49 discussed entropy optimization rate in MHD nanofluid flow through micro parallel networks.Khan et al. 50scrutinized irreversibility in rotational viscous nanofluid flow with mixed convection and radiative heat flux.Huminic and Huminic 51 reported irreversibility in hybrid nanomaterials flow.
The literature cited in above transpires that irreversibility in Darcy Forchheimer flow of Sutterby nanomaterial due to stretched cylinder with porous walls in existence of Arrhenius kinetics, Joule heating and chemical reaction is not examined till now.In order to fill this gape, motivation here is to investigate the irreversibility in radiative Sutterby nanoliquid flow by stretchable cylinder.Heat transport characteristics are scrutinized through Joule hating and radiation effects.Furthermore, Brownian movement and thermophoresis diffusion impacts are accounted.Mass transfer characteristics are reported through activation energy.Utilizing thermodynamics second law physical description of irreversibility is analyzed.

Problem description
In current inspection, incompressible steady flow of Sutterby nanomaterial by stretchable cylinder is considered.The flow is taken along axial z-direction and radial direction is taken normal to z-direction.The cylinder stretches with velocity U w = U 0 z l in the axial direction due which flow generates.Magnetic field of constant strength B 0 is imposed vertical to the flow.Induced magnetic field effect is neglected for small Reynolds number.Effects of Arrhenius kinetics, Joule heating and radiation have been incorporated in thermal and mass concentration equations.Boundary layer conventions are accounted in development of flow governing model equations.Schematic flow diagram with boundary restrictions is depicted in Fig. 1.The governing equations signifying the flow under above norms are as follows 52,53 ; with (1) ∂u ∂r

Entropy modeling
It is defined as 52,53 ; non-dimensional form of entropy generation rate is Bejan number is here the characteristic entropy rate,

Numerical solution and discussion
Here NDsolve code in Mathematica is executed to solve the dimensionless system of equations.Impact of sundry variables on Sutterby nanofluid velocity (f ′ (η)) , thermal field (θ(η)) , mass concentration (φ(η)) EG rate (S G ) and Bejan number (Be) are scrutinized by plotting.Engineering quantities are evaluated numerically.Table 1 is labeled to ensure the correctness of present numerical approach.This table demonstrated the comparison of gradients of temperature θ ′′ (0) versus different Pr values while influence of remaining variables is neglected.Clearly the results are in good agreement.

Velocity
Figure 2 delineates Sutterby fluid parameter (β 1 ) impact on f ′ (η) .For raising values of β 1 velocity decreases.Physically for raising β 1 relaxation time increases as a result viscous effects dominants therefore additional resistance is offered to fluid particles thus f ′ (η) decreases.Figure 3 discovers influence of Forchheimer variable

∂T ∂r r=R
, q m = −D B ∂C ∂r r=R . ( Cf z (Re z ) (Fr) on fluid velocity.Physically, for hiking Fr estimations, f ′ (η) decreases.For larger values of Fr drag surface force increases, thus f ′ (η) decreases.Figure 4 shows influence of γ 3 on f ′ (η).here, material velocity upsurges via rising γ 3 .For higher γ 3 cylinder radius decreases as a result fluid contact area between fluid particles and surface of cylinder decreases due to which less resistive force of surface is offered to fluid particles so velocity upsurges.Figure 5 describes influence of surface porosity on f ′ (η) .Velocity decreases for improvement in .For increasing values of kinematic viscosity of porous medium upsurges, therefore f ′ (η) curves decays.Behavior of f ′ (η) versus magnetic variable is captured in Fig. 6.Velocity diminishes versus higher M .Since through higher M Lorentz force enhances which acts in the opposite direction of flow and thus velocity decays.From Fig.

Temperature
Figure 8 gives inspiration of Ec on θ(η).For enlargement in Ec temperature upturns.Since for higher Ec , internal energy of fluid boosts, consequently kinetic energy of system increases as a result inside friction of tiny solid particles enhances and extra heat supplied to the system, resultantly θ(η) increases.Variation in thermal fluid versus γ 3 is captured in Fig. 9. Clearly θ(η) increases for raising values of γ 3 .Figure 10 shows the variation in   θ(η) against M .For greater estimations of M temperature rises.In fact more heat is added in the system when magnetic variable amplifies due to resistive force and so θ(η) increases.Influence of Nb on θ(η) is checked in Fig. 11.Temperature improves through higher Nb.Random movement of tiny particles increases with in fluid against higher Nb , consequently inter collision of nanoparticles increases and thus extra heat generates, therefore θ(η) improves.Figure 12 is designed to check the outcome of Nt on θ(η).For raising Nt thermal field boosts.Physically, when thermophoresis force increases rate of shifting of solid tiny particles from hot to cold region improves, so θ(η) escalates.Figure 13 articulates the influence of Pr on thermal field.It is perceived here that via higher Pr thermal field declines.Physically through rising Pr fluid thermal diffusivity diminished and thus θ(η)

Concentration
Figure 15 is drafted to show influence of δ 1 on φ(η) .It is observed that for an escalation in δ 1 values φ(η) upsurges.Figure 16 depicts E 1 effect on φ(η) Concentration upsurges for raising values of E 1 .For higher E 1 Arrhenius func- tion decreases thus concentration increases.Figure 17 portrays effect of γ on φ(η) .Concentration decreases for raising γ .For higher γ destructive rate of reaction increases and thus liquid species liquefy more successfully, so   φ(η) decreases.Figure 18 gives impact of γ 3 on φ(η) .For enlargement in γ 3 concentration increases.Figure 19 expresses the Nb influence on φ(η) .Nanofluid mass concentration declines through rising Nb .Physically, inter- molecular collision boosts against higher Nb due to arbitrary movement.Resultantly more heat generates thus temperature increases and φ(η) decreases.www.nature.com/scientificreports/transfer rate is examined in Table 4. From this table it is observed that Sherwood number improves for higher δ 1 , γ , n, Nb, Nt and Sc .Opposite effect is seen for E 1 and γ 3 .

Conclusions
The prime objective of this article is to examine the irreversibility in steady magnetized flow of Sutterby nanofluid caused by stretched cylinder.The features of Arrhenius kinetics, Joule heating, chemical reaction, Darcy Forchheimer, surface permeability and thermal radiation are accounted in development of mathematical governing   www.nature.com/scientificreports/equations.Through thermodynamics 2 nd law total irreversibility is modeled.Numerical and graphical solutions are constructed through RKF-45 in Mathematica package.Main findings are itemized as: • Velocity decreases for higher porosity parameter and magnetic variable.
• Temperature have opposite behavior for magnetic variable and Prandtl number.
• For higher thermophoresis and Brownian movement variable temperature enhances.www.nature.com/scientificreports/ • Concentration decreases for Brownian movement variable while it increases for thermophoresis variable.
• For higher Schmidt number concentration decreases.
• EG increases for higher Brinkman number, temperature difference ratio variable and magnetic variable.
• Bejan number decays versus higher Brinkman number and magnetic variable while it enhances for tempera- ture difference ratio variable.

Figure 1 .
Figure 1.Coordinate system and flow diagram.
Figure 20 explores the behavior of φ(η) versus Nt .Fluid concentra- tion increases for higher values of Nt .Physically, thermophoretic force which shifts fluid particles from hot to cold region increases for enlargement in Nt thus φ(η) increases.
Figure 21 is plotted to study influence of Sc on φ(η) Here, φ(η) decreases via increasing values of Sc.Since mass diffusivity reduces for larger Schmidt number and thus φ(η) is diminished.

Figure 26 .
Figure 26.S G via M.

Table 1 .
Comparison of results for gradient of temperature.